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Quasiparticle Excitations outside the Vortex Cores in MgB2 Probed by Muon Spin 

Rotation 
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The magnetic penetration depth A in the mixed state of MgB2 has been determined micro- 
scopically from the field distribution around magnetic vortices probed by muon spin rotation. 
Both the temperature and magnetic field dependence of A strongly suggest the presence of 
delocalized quasiparticles around the vortex cores. In particular, the effect of Doppler shift has 
been clearly observed as a finite gradient of the field dependence of A, strongly suggesting an 
anisotropic order parameter with the region of a vanishingly small energy gap. 
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The revelation of superconductivity in a binary inter- 
metallic compound, MgB2, has attracted much interest 
because it exhibits a transition temperature (~ 39 K) 
almost two times higher than those of all other inter- 
metallic superconductors known to date.^^ The most in- 
teresting issue associated with this compound is whether 
or not it belongs to the class of conventional Bardeen- 
Cooper-Schrieffer (BCS) type superconductors. To date, 
most experimental results have favored a phonon me- 
diated superconductivity. The boron isotope effect,^-' 
photoemission spectroscopy,'^-' ^^B-NMR,'') Raman spec- 
troscopy,^^ tunneling measurements, ^"^^ and optical con- 
ductivity data^*^-* are all consistent with the conventional 
BCS s-wave pairing. On the other hand, calculations of 
the band structure and the phonon spectrum predict a 
double energy gap,^^'-^^) a larger gap attributed to two- 
dimensional Px-y orbitals, and a smaller gap attributed 
to three-dimensional Pz bonding and antibonding or- 
bitals. Experimental results of specific heat measure- 
ments,^^' point-contact spectroscopy, photoemis- 
sion spectroscopy,^^-' scanning tunneling spectroscopy^^-' 
and penetration depth measurements^^-* have supported 
this scenario. 

It must be noted that earlier experiments includ- 
ing muon spin rotation (/iSR)^^-* and ac susceptibil- 
j^yig, 20) performed on polycrystalline samples revealed a 
quadratic behavior of A at low temperatures, from which 
they inferred the presence of line nodes in the order pa- 
rameter. However, the recent /xSR analysis has demon- 
strated that such behavior can also be explained by as- 
suming the double energy gap (without resorting to the 
line nodes). Thus, the temperature dependence of A 
provides limited information for determining the struc- 
ture of the order parameter. This situation can be im- 
proved by studying the magnetic field dependence of A, 
where the quasiparticle excitation is controlled by the 
Doppler shift^^^ which is independent of the thermal ex- 
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citation. In order to obtain more detailed information 
on the order parameter in MgB2, we have observed the 
field dependence of A over a wide range of magnetic field 
up to 5 T. 

The magnetic penetration depth A is determined by 
the quasiparticle excitations outside the vortex cores, 
and thus it provides an excellent measure of the struc- 
ture of the order parameter in the fiux-line lattice (FLL) 
state. The /LtSR technique is a powerful microscopic tool 
for obtaining the fundamental length scale such as A in 
the bulk type II superconductors. Implanted muons ran- 
domly probe the local magnetic fields induced by the 
FLL, yielding the spectral density n{B) which is directly 
related to the spatial field distribution B{r). While the 
/iSR spectra in a single crystalline specimen can be com- 
pared directly with n{B) calculated from the spatial field 
distribution B{r), the spectra in a polycrystalline speci- 
men are subject to the modulation of line shape due to 
various kinds of inhomogeneity. Even in such a situa- 
tion, the time-dependent /zSR spectra is approximately 
described by a Gaussian damping exp(— tT^t^/2) with a 
being primarily determined by the second moment of 
n{B). 

In this letter, we report on the temperature and field 
dependence of A deduced from those of a in polycrys- 
talline MgB2. A special precaution has been taken 
in selecting the field for the temperature scan, in or- 
der to avoid the effect of random flux pinning near the 
lower critical field Hd- Analysis based on the two-gap 
model yields Ai = 4.9(1) meV and A2 = 1.2(3) meV. 
The penetration depth at T = K is estimated to be 
103.9(1.0) nm. We also found that A at T ~ 10 K exhibits 
a significant increase with almost linear dependence on 
the applied magnetic field, which can be understood by 
considering the Doppler shift of the quasiparticle excita- 
tion associated with the anisotropic order parameter. 
However, the gradient against the field is considerably 
small compared with that in c?-wave superconductors. 
These results indicate the existence of excess quasiparti- 
cle excitations outside the vortex cores in MgB2, strongly 
suggesting that there is an anisotropic structure in the 
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order parameter with a nodal region smaller than that 
for the d-wave pairing. 

The polycrystalline sample of MgB2 used in this ex- 
periment had a surface area of ~ 50 mm^. The super- 
conducting transition temperature Tc determined from 
resistivity and susceptibility measurements was 38.5 K.-^^ 
nSK experiment was performed on the M15 beamline at 
TRIUMF which provides a muon beam with the mo- 
mentum of 29 MeV/c. A muon- veto counter system was 
adopted to eliminate positron events from muons which 
missed the sample so that the relative yield of such events 
was less than 5% of the total positron events. An exper- 
imental setup with a high time resolution was employed 
to measure the transverse field (TF-) fiSR time spectra 
up to 5 T. The sample was field cooled at the measured 
magnetic fields in order to eliminate the effect of flux 
pinning. 

Since the muons stop randomly on the length scale of 
the FLL, the muon spin precession signal P{t) provides a 
random sampling of the internal field distribution B{r), 

/oo 
n{B)exp{ij^Bt)dB, (1) 
-OO 



n{B) = 



dr 
dB' 



(2) 



where 7^ is the muon gyromagnetic ratio (= 
2nx 135.53 MHz/T), and n{B) is the spectral density de- 
termined by the local field distribution. These equations 
indicate that the real amplitude of the Fourier trans- 
formed muon precession signal corresponds to the spec- 
tral density n{B). The London penetration depth in the 
FLL state is related to the second moment ((Ai?)^) of 
n{B) reflected in the /xSR line shape. For polycrys- 
talline samples, a Gaussian distribution of local fields is 
a good approximation, where 



P{t) ~ exp(-(7^r/2)exp(i7^iJf) 



J.y^iiABT), 



(3) 

(4) 



with H being the external magnetic field. For the case 
of an ideal triangular FLL with the isotropic effective 
carrier mass m* and a cutoff K « 1.4/^„ provided by 
the numerical solution of the Ginsburg-Landau theory, 
the London penetration depth A can be deduced from a 
using the following equation:^^' 

a [pis-^] = 4.83 X 10^(1 -h)[l + 3.9(1 - hfj^/^X'^ [nm], 

(5) 

where h = H/Hc2- It should be noted that eq. (5) pro- 
vides the field dependence of a when A is a constant. 
Here, A is related to the superconducting carrier density 
rig as follows: 



A^ = 



(6) 



indicating that A is enhanced upon the reduction of rig 
due to the quasiparticle excitations. 

Figure 1 shows the fast Fourier transforms (FFT) of 
the muon precession signal in MgB2 for different temper- 
atures at if ~ 0.5 T, where the real amplitude of FFT 
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Fig. 1. FFT spectra of ^SR signal in MgB2 at ~ 0.5 T under 
several temperatures. 



corresponds to n{lS) in the FLL state convoluted with 
additional damping due to small nuclear dipolar fields. 
As is most explicit in the FFT spectrum at T = 8.5 K, 
the line shape is characterized by a broad peak near 
B — ~ and a satellite at lower field (with a small 
background peak right at i3 = iJ). This is a typical fea- 
ture often seen in polycrystalline powder samples with 
a lower Meissner fraction, where the sample consists of 
both superconducting and normal domains with the typ- 
ical size of a few microns (> A). The magnetic field 
distribution at the normal domain is shifted to a high 
field because of the demagnetization associated with the 
Meissner effect in superconducting domains. Thus, the 
peak at the lower field corresponds to the signal from the 
superconducting domains. 

Considering the double peak structure in the FFT 
spectra in Fig. 1, we adopted two components with the 
following empirical form for analyzing the data in the 
time domain: 



A^P{t) = ^ Aj(t)exp[i(7^Bjt- 
A,-(t)=A,exp[-((7,t)V2], 



(7) 



(8) 



where ^0 is the total positron decay asymmetry, Aj {j = 
1,2) the partial asymmetry for respective components, 
Bj the central frequencies, (j) the initial phase, and aj 
the muon depolarization rates. 
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Fig. 2. Temperature dependence of the muon depolarization rate 
(Ti under H = 0.5 T. The sohd curve is the result of fitting by 
eq. (9). 



We show the temperature dependence of the muon de- 
polarization rate ai (which corresponds to the linewidth 
of the lower frequency peak in Fig. 1) in Fig. 2. These 
data were obtained under a field ~ 0.5 T which is well 
above Hd ~ 10~^ T. Following the results of specific 
heat,"'^^' "'^^^ Raman^^ and /iSR measurements,^ the data 
were fitted by the two-gap model^^' which is described 
as follows: 

a{T) = a{0) - w ■ 6cr{Ai,T) - {1 - w) ■ Sa{A2,T), (9) 
2a(0) 



= /(£,r).[l-/(£,r)]&, 



(10) 



where w is the ratio of gap energy between the two gaps, 
ks the Boltzmann constant, f{s,T) the Fermi distribu- 
tion of quasiparticlcs. and A(T) the BCS gap energy.^"''' 
The solid line in Fig. 2 is the best fit result with Ai = 
4.9(1) meV, A2 = 1.2(3) meV, w = 0.84(3) and Tc = 
36.48(3) K. The dotted line shows the result of fitting 
using the values of Ai, A2 and w in ref. 21. Although 
our result shows reasonable agreement with the earlier 
one,^^-* the value of A2 is considerably smaller than the 
reported value of 2.6(2) meV. 

Figure 3(a) shows the magnetic field dependence of the 
muon spin relaxation rate cri with an inset showing the 
low field region {H < 0.6 T). As seen in the earlier re- 
port, the effect of random flux pinning is observed as a 
peak of (71 in the low field region (see inset of Fig. 3(a)). 
However, ui decreases with increasing field above H = 
0.1 T, indicating that the distortion of FLL is reduced 
by increasing inter- vortex interaction and that the depo- 
larization is predominantly determined by the intrinsic 
n{B). The dashed line in Fig. 3(a) shows the fitting re- 
sult by eq. (5) with Hc2 = 12.5 T as determined by resis- 
tivity measurements.^®) (Here, we discuss the magnetic 
field dependence of ui for the data above = 0.5 T to 
avoid the remnant effect of flux pinning at lower flelds.) 
Compared with the field dependence of ai, the dashed 
line does not reproduce the data, indicating that A is not 
a constant but it increases with increasing external field. 
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Fig. 3. Magnetic field dependence of (a) the muon depolarization 
rate g\ and (b) the magnetic penetration depth A at T ~ 10 K. 
Dashed and solid lines are the fittings using (a) eq. (5) with a 
constant A or A proportional to A(0)[1 -|- r; • ft] and (b) eq. (12). 
H^2 at this temperature is approximately 12.5 T. 



The field dependence of A estimated from eq. (5) is 
shown in Fig. 3(b). It clearly exhibits a strong field de- 
pendence where A(/i) increases almost linearly with h. 
This is similar to the cases of YNi2B2C,2'^) NbSe2^*) and 
high- Tc cupratc superconductors, where the increase 
of A is attributed to the anisotropic order parameters and 
the associated nonlinear effect due to the Doppler shift of 
the quasiparticlcs in the nodal region (A(fc) ~ 0).^^^ The 
field dependence of A is expected to be stronger when the 
phase space satisfying A(fc) ~ has larger volume. A 
fltting by the relation 



A(/7) = A(0)[l+r/-/i] 



(12) 



provides a dimensionless parameter r\ which represents 
the strength of the pair breaking efli'ect. We obtain = 
1.27(29) with A(0) = 116.1(6.2) nm which is shown as the 
solid line in Fig. 3(b). The obtained value of r\ is inter- 
mediate between that in YNi2B2C and NbSe2 (e.g., 77 = 
0.97 at 0.2rc,2^) 7] = 1.61 at 0.33Tc,2«) respectively) and 
is smaller than those in d-wave superconductors (e.g., r\ 
= 5.5 ~ 6.6 for cuprates^^^). The solid line in Fig. 3(a) 
is the result of fltting with the relation of eq. (5), with A 
represented by eq. (12). 

Our result on the temperature dependence of A is 
qualitatively consistent with earlier results, ^-'^^ suggest- 
ing that the order parameter in MgB2 may be effectively 
described by adopting the two-gap model. However, the 
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observed field dependence of A is not expected for the 
isotropic order parameter irrespective of the muhiplic- 
ity of the band structure and the associated gap energy. 
Considering that the current rcsuh on the field depen- 
dence of A was obtained at T ~ 10 K, this energy scale 
of e = fcsT ~ 1 meV places an upper boundary on 
the smaller gap energy A2 in order to explain the ob- 
served effect of the Doppler shift. Since our estimation 
for A2 = 1-2(3) meV is very close to e, the observed 
if-linear behavior of A may be attributed to the quasi- 
particle excitations in the vicinity of the smaller gap. 
While this cannot be distinguished from the case of a 
nodal structure in the order parameter (i.e., considering 
a region where A(A;) <C e), our result is clearly incon- 
sistent with the two-gap model with A2 » e. On the 
other hand, it is also quite unlikely that the d-wave pair- 
ing is realized in MgB2, because the coefficient rj is much 
smaller than those in high- Tc cuprates. The recent obser- 
vation that the order parameter in YNi2B2C (where the 
pairing symmetry has been identified as s-wave^^)) has 
point nodes^^^ exhibits a good correspondence with the 
intermediate value of ?7(~ 1) obtained by /uSR,^'') sug- 
gesting that there is a similar situation in MgB2. Thus, 
the present ^SR result leads us to conclude that the or- 
der parameter in MgB2 has a structure with an energy 
gap smaller than e ~ 1 meV. The field dependence of 
A measured at a much lower temperature would provide 
more useful information to distinguish the anisotropic 
order parameter from the isotropic one described by the 
two-gap model. 

In summary, we have performed TF-/xSR measure- 
ments in MgB2 to obtain the temperature and magnetic 
field dependence of the penetration depth A and the asso- 
ciated spin relaxation rate Our result is perfectly in 
line with the presence of an anisotropic order parameter 
with a nodal structure, and it sets an upper boundary 
e 1 meV for the smaller gap energy in the two-gap 
model. The magnetic field dependence of A exhibits a 
linear dependence on the external field up to 5 T with 
the gradient r] being considerably smaller than that in 
d-wave superconductors, which may disfavor the occur- 
rence of rf-wave pairing in MgB2. 
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